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Abstract 

We study the temporal aspects of laser-assisted extreme ultraviolet (XUV) photoionization using attosecond pulses 
of harmonic radiation. The aim of this paper is to establish the general form of the phase of the relevant transition 
amplitudes and to make the connection with the time-delays that have been recently measured in experiments. We 
find that the overall phase contains two distinct types of contributions: one is expressed in terms of the phase-shifts of 
the photoelectron continuum wavefunction while the other is linked to continuum-continuum transitions induced by the 
infrared (IR) laser probe. Our formalism applies to both kinds of measurements reported so far, namely the ones using 
attosecond pulse trains of XUV harmonics and the others based on the use of isolated attosecond pulses (streaking). 
The connection between the phases and the time-delays is established with the help of finite difference approximations 
to the energy derivatives of the phases. The observed time-delay is a sum of two components: a one-photon Wigner-like 
delay and an universal delay that originates from the probing process itself. 
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The dynamics of photoionization can now be explored 
with unprecedented time resolution thanks to high-order 
harmonic-based sources that deliver pulses of XUV radia- 
tion with duration in the attosecond range. Recent mea- 
surements performed with single attosecond pulses have 
shown the existence of an unexpected time-delay between 
the single-photon ionization from the 2s and the 2p sub- 
shells of Ne atoms in gas phase Jl, . The "streak camera" 
technique used in these experiments [2] implied nontriv- 
ial ejection times of the photoelectrons, depending on the 
sub-shell from which they originate. Similar delays be- 
tween the ejection times from the 3s and 3p sub-shells 
in Ar have been measured also using trains of attosecond 
pulses [3] , with the help of another technique based on in- 
terferometry called RABBIT (Reconstruction of Attosec- 
ond Beating By Interference of Two-photon transitions) 
[3HS]- In both cases, delays of several tens of attoseconds 
have been measured. As photoionization is one of the most 
fundamental processes in light-matter interactions, these 
results have motivated a large number of theoretical inves- 
tigations [THE]- 

The two kinds of measurements share many similarities 
since they involve a laser-assisted single-photon ionization 
process and they rely on a phase-locked IR laser field to 
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probe the temporal aspects of the XUV photoionization. 
However, they differ in the analysis used to determine the 
time-delays and in the range of IR laser intensity. 

The motivation of the present paper is to present an 
unified theoretical analysis of these processes. To achieve 
this goal, we shall expose first the theoretical background 
which has conducted us to conclude in [5, that in inter- 
ferometric measurements, the measured delays arise from 
the combination of two distinct contributions: One is re- 
lated to the electronic structure of the atomic target while 
the other is induced by the measurement process itself. 
The first one can be identified as a "Wigner time-delay" 
[m [15] that is directly related to the energy dependence 
of the different phase-shifts experienced by the photoelec- 
trons ionized from distinct sub-shells in atoms. The other 
contribution is induced by the IR laser field that is used 
to probe the photoionization process. This latter contri- 
bution results from the continuum-continuum transitions 
induced by the probe IR laser field in the presence of the 
Coulomb potential of the ionic core. When simplifying 
the analysis to the cases when the process is dominated 
by the asymptotic form of the relevant second-order ma- 
trix elements, a characteristic measurement-induced delay 
can be identified, that is independent from the details of 
the electronic structure of the ionic core. This shows how 
the experimental signal can be related to the temporal dy- 
namics of one-photon ionization. 

Regarding the streaking measurements realized with a 
single attosecond pulse of XUV radiation |1], the experi- 
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Figure 1: (a) Laser-assisted photoionization by an attosecond pulse train, corresponding to odd XUV harmonics, H2q—i, where q is a positive 
integer. The sideband S2q can be reached by either absorbing H'zq—i and then absorbing a IR laser photon uj or by absorbing H29+1 and 
then emitting ui. The sideband signal, S2q, oscillates as a function of the subcycle-delay, r, between the attosecond pulses and the IR laser 
probe field, (b) Laser-assisted photoionization by a single attosecond pulse, corresponding to a broad XUV continuum. To first order, the 
electron is ionized by absorbing one XUV photon, Q, resulting in a wave packet centered at e = + f2o- To second order, the electron 
may absorb an additional laser photon uj resulting in an upshifted wave packet centered at e = -|- Qq + or it may emit a laser photon 
resulting in a downshifted wave packet centered at e = -|-f!o —(^- The interference of these three wave packets leads to a modulation, Afc, of 
central momentum of the photoelectron as function of the subcycle-delay, r, between the laser field and the attosecond pulse, (c) Two-photon 
XUV-IR Above-Threshold Ionization from an initial bound state with energy e^, to a final state with energy e. 



mental data were obtained for IR field intensities signifi- 
cantly higher than those obtained with attosecond pulse 
trains [5]. Understandably, the questions related to the 
role of the probe IR field on the photoelectron dynamics 
in streaking measurements have motivated several theo- 
retical studies [7HI3]: see also the earlier papers: [T6HT9] . 
Then, a natural issue arises which is to determine to what 
extent the "streaking delays" so obtained differ from those 
derived from the interferometric data. Although both the 
experimental techniques and the theory treatments differ, 
it is of interest to compare the two approaches. Indeed, 
as we shall show below, a link can be found when re- 
ducing the laser intensity of the streaking field so that 
one reaches the domain of applicability of the recently 
developed Phase-Retrieval by Omega Oscillation Filter- 
ing (PROOF) scheme, [20]. An interesting outcome of 
our analysis is to show the importance of the long-range 
Coulomb potential for understanding the absolute time- 
delays in the streaking experiments as well. 

The interpretation of the attosecond delays in pho- 
toionization relies on our ability to determine the phases 
of the relevant transition amplitudes. Thus, before go- 
ing into the details of the derivation of such phases, we 
shall outline the main features of the two techniques in 
Section [2] Then, Section [3] is devoted to the presenta- 
tion of the general expressions for two-color, two-photon, 
complex transition amplitudes that are relevant for Above- 
Threshold Ionization (ATI) in single-active electron sys- 
tems. The theoretical background is based on a perturba- 
tive approach and the emphasis will be on the derivation of 



a closed-form approximate expression that is of interest for 
evaluating the phase of the amplitudes. The basis of exact 
computations in hydrogen will be outlined, and a simpli- 
fied classical treatment will be presented. Applications to 
the determination of the relation between the phases and 
the time-delays is presented in Section|4j Here we consider 
first ionization by an attosecond pulse train and then by a 
single attosecond pulse, in the presence of a relatively weak 
IR field. This discussion provides an interesting connec- 
tion between the two types of measurements. Section [5] 
contains a comparison of the results extracted from the 
approximate evaluation of the delays to the ones deduced 
from exact calculations performed in hydrogen from dif- 
ferent initial states. Also, we present our conclusions and 
perspectives. 

2. Laser- Assisted XUV Photoionization: Attosec- 
ond Pulse Train vs. Single Attosecond Pulse 

The principle of the measurements of the delays using 
an attosecond pulse train is illustrated in Fig. [l](a), which 
represents schematically the ionization of an atom in the 
simultaneous presence of a set of several XUV (odd) har- 
monics and of the IR field, used to generate the harmonics 
(atomic units will be used throughout the paper, unless 
otherwise stated). In the time domain, both pulses are 
"long", i.e. the IR laser pulse is multi-cycle, with typical 
duration of a few tens of femtoseconds, and the XUV har- 
monic field is constituted of a train of attosecond pulses 
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(or equivalently of a comb of coherent odd harmonic fre- 
quencies (2(7 + H2q+i)- Under these conditions, the 
photoelectron spectrum consists of equidistant lines sepa- 
rated by 2uj that are associated to one-photon ionization 
of the target by each harmonic. In-between these lines are 
sidebands associated to two-photon transitions involving 
the absorption of one harmonic and the exchange of one 
IR photon. The signal intensities, S2q^ of the sidebands la- 
belled 2q vary periodically with the delay t between the IR 
and the harmonic pulses, according to a generic expression 
that involves the phases of the fields together with atom- 
dependent contributions: 
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a -I- /3C0S[2WT - A02g - ^02q], 



(1) 



where A(/)2g = {(j)2q+i — 4>2q-i) is the phase difference be- 
tween the consecutive harmonics -ff2g+i and i?2g-i and 
A92q is an intrinsic atomic quantity, associated to the dif- 
ference of the phases of the transition amplitudes associ- 
ated to the distinct quantum paths leading to the side- 
band [5]. 

To make clearer the connection between the above phase 
differences and the time-delays we shall discuss here, it is 
convenient to rewrite the formula in Eq. ([!]) under the 
form: 
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a + (3 cos[2a;(r - T2q - Tg)] 



(2) 



where T2q = A02g/2w is a finite difference approximation 
to the group delay GD = dcj)/dft of the harmonic radia- 
tion at the considered frequency, Q w 2quj, as presented 
in refs. [31] [55]. On the other hand, Tg = A92q/2uj is an 
intrinsic time-delay associated to the atomic phase differ- 
ence, A02g- As reported in [3] and as we shall describe 
in more details here, the determination of Tg gives access 
to the temporal dynamics of atomic photoionization. Be- 
fore closing this brief presentation of the RABBIT scheme, 
we stress that the intensities of both fields must be kept 
moderate, so that the phases of the transition amplitudes 
associated to the sidebands can be derived from a standard 
time- dependent perturbation theory calculation, limited to 
second-order. 

As represented schematically in Fig. [l] (b), streaking 
relies on the ionization of the atom by a single attosec- 
ond pulse, in the presence of a few-cycle IR pulse. One 
requirement to realize streaking is that the effective dura- 
tion of the attosecond pulse has to be significantly shorter 
than the IR pulse cycle [5] (more rigorously it is actually 
the spectral bandwidth of the attosecond pulse that must 
be larger than the probe photon frequency). The mea- 
surement consists then in recording the momentum, kf{t), 
of the ejected photoelectron, as deflected by the instanta- 
neous IR probe vector potential, Ai^{t), so that its wave 
vector is given approximately by: 



kf Kik-A^{t), 



(3) 



where k is the field-free momentum. We note that this 
simple relation is being derived by assuming that the pho- 
toelectron does not experience the effects of the residual 



Coulomb potential of the ionic core [5|. In this article, we 
will recover this streaking phenomenon using an interfero- 
metric interpretation, which includes the full effect of the 
ionic core, thereby, obtaining the correct absolute delay of 
the momentum modulation relative to the probe field. 
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Figure 2: Sketch of an attosecond photoionization time-delay ex- 
periment between two different initial states with energies, and ej . 
The attosecond XUV field, O, ionizes the electron in the presence 
of a phase-locked IR laser field, u!. Since the same single attosecond 
pulse (or attosecond pulse train) is used to promote the electrons 
from either state, the observed delays of the modulation of the cen- 
tral momenta (or sidebands), ATij{Q), are directly related to the 

atomic delay difference T^''(f2) — r^^''(0) defined in Eq. (j^Jl, which, 
in turn, is related to the difference of the corresponding two-photon 
matrix element phases. 

If the XUV field frequency is high enough to ionize 
electrons from either one of two valence sub-shells of an 
atom (typically the 2s and 2p states in Ne atoms) one 
can record two distinct streaking traces corresponding to 
the two photoelectron lines associated to each sub-shell. 
This situation is schematically displayed in Fig. 2. A de- 
lay between the ejection times of the photoelectrons from 
the different atomic sub-shells can be determined by com- 
paring the corresponding streaking traces, corrected from 
possible biases introduced by the experimental procedure 
[1]. The same idea can be applied to attosecond pulse 
trains, by using the interferometric set-up [3]. We turn 
now to the presentation of the theoretical background that 
is common to the two kinds of techniques, in the limit of 
weak IR probe fields. 
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3. Phases of Laser- Assisted Photoionization Tran- 
sition Amplitudes 

3.1. Two-photon Above-Threshold Ionization 

A representative laser-assisted photoionization transi- 
tion is depicted in Fig. [l] (c) : It displays the sequential ab- 
sorption of one XUV harmonic photon with frequency fi, 
followed by the absorption of one IR laser photon with fre- 
quency bj. It corresponds to the lowest-order perturbative 
amplitude for an Above- Threshold Ionization (ATI) pro- 
cess observed when the XUV frequency is larger than the 
ionization energy of the system: D, > Ip. Obviously, other 
quantum paths are contributing to this type of two-color 
ionization process, e.g. the IR photon can be absorbed be- 
fore the XUV photon, but ATI amplitudes of the former 
type are dominant in the class of experiments considered 
here. 

For two fields with the same linear polarization e, it is 
natural to choose this direction as the quantization axis z, 
and the matrix element associated to the path shown in 
Fig. [l] (c), is of the form: 



M{k; + O) ^ -^E^En lim 



im V 



k \ z \ V \ z \ i] 



(4) 



where Eq and E,^ are the complex amplitudes of the har- 
monic and IR laser fields, respectively; ipni,ei,mii'r) = {r\i) 
is the initial state wavefunction, with negative energy 
and = {'r\k) is the final state wavefunction with pos- 

itive energy = = ei -\- Q. -\- ui . The sum over the 

index v runs over the whole spectrum (discrete plus con- 
tinuous) of the atom. The partial wave expansion of the 
final state wavefunction is: 

^-(rO = (8^)=^/^ *^e-''-('=)r£,MWi^L,M(r)i?fe.L(r), 

where the Y^^m are spherical harmonics, the Rk,L(j') are 
(real) radial wavefunctions normalized on the energy scale 
and rjL{k) are the phase-shifts. We note that this wave- 
function behaves asymptotically as the superposition of a 
plane wave plus an ingoing spherical wave, as required to 
treat photoionization Thus, the phase-shifts 7]L{k) 

account for the phase difference between the free motion 
of a plane wave and that of a photoelectron wave ejected 
from an atomic bound state. 

The angular dependence of the matrix element can be 
factorized out for an initial state 



and with z = •\/47r/3 r Yi Q{f), it becomes: 
M{k, + f}) - — (87r)3/2 ^^-^^ 
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L,M 

X 5](rL,M|>l,o|n,^)(>A,Mln,o|>£.,™.) 

X TL,xA{k,e, + n), (5) 

where the angular momentum components of the interme- 
diate states are labelled (A, /i) and the quantity, Tl.a,^; [k; e^- 
Q), is the radial part of the amplitude. The span of ac- 
cessible angular momentum states in the intermediate and 
final states is governed by the dipole selection rules: One 
has X = iizLl; L = £i,£i±2 and M — n = rui respectively. 
The explicit form of the radial amplitude T^^x^i. (fc; -I- Q) 
is: 



{Rk,L\r\R,,x){R,,xMRn,,e,) 



lim de (^'^^L\ r\R^',x){R.',x\r\Rn,A) 

lo 



(6) 



where we have separated the contributions of the discrete 
and continuous spectra. 

Since the frequency of the XUV harmonic is larger than 
the ionization potential of the atom, f2 > /p = |ei|, it is also 
larger than the excitation energies of the atom, H. > e^ — ei. 
Accordingly, the denominators of the terms in the sum 
over the discrete states are positive and relatively large, 
which makes the overall contribution of these terms sig- 
nificantly smaller than that of the continuous spectrum. 
In the integral running on the continuous spectrum of en- 
ergies — the denominator becomes zero at the 
energy — = -|- fi. Taking the limit e ^ 0^, the 
integral becomes: 

lim de , ^^k^L\r\R.',x){R.',xV\Rna:) 

e-)-o+ Jo " ei + VL- -I- ie 

f +00 



= v 



{RKLV\R^'.x){Rn',x\r\Rn^.u) 

-i^{Rk,L\r\R.,x){R.,x\r\Rn^A)- 



(7) 



where the first term is the Cauchy principal value of the 
integral, which turns out to be real, and the second term 
is purely imaginary. The latter is associated with a two- 
step transition as it contains the product of the one-photon 
ionization amplitude towards the state of energy = ti-\- 
r2, times the continuum-continuum transition amplitude 
from e„ towards the final state of energy — e„ + w 

that is reached upon the absorption of the IR photon lo. 

The overall phase of the radial matrix element, Eq. ([6]), 
is thus governed by the ratio of the imaginary term in 
Eq. ([?]) to the sum of the integral principal part in the same 



equation plus the contribution of tlie discrete spectrum 
contained in Eq. (|6| . Accurate computations of such am- 
phtudes and phases represent a formidable task for most 
atomic systems. This entails to rely on approximate rep- 
resentations of the atomic potential for each angular mo- 
mentum dependent state [551 [IS] or to use many-electron 
techniques [571 HH] • There is however the notable excep- 
tion of hydrogenic systems, where "exact" calculations of 
these amplitudes are feasible [5^32) . see below. It is thus 
of importance to derive an approximate treatment which 
should allow to get correct estimates of the phases of in- 
terest to address the questions of the time-delays. 

3.2. Asymptotic approximation for 2-photon ATI matrix 
elements 

Let us re-express the radial amplitude in terms of the 
first-order perturbed wavefunction denoted Pn^xir): 



(8) 



The function Pk,x(j) solves the inhomogeneous differential 
equation: 



(9) 



where H\ is the radial atomic Haniiltonian for angular mo- 
mentum A. The fully developed eigenfunction expansion 
of PK,\{r) can be identified using Eqs. Q and Q: 
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(10) 



We note that it describes the radial part of the intermedi- 
ate photoelectron wave packet created upon absorption of 
the XUV photon f2, before absorbing the IR laser photon 
(jj. 

The essence of the approximate treatment that we have 
implemented to get an estimate of TL^x,i.i(k;ei -I- ri), is 
based on using the asymptotic forms of both the final 
state function Rk^hif) and of the perturbed wavefunction 
PK,\{f) for large values of their radial coordinate. This 
is a priori justified by the fact that we are interested in 
the phases of the amplitudes which are governed by the 
asymptotic behavior of these functions. As an additional 
verification, we will compare the predictions of the model 
to those derived from an exact treatment in Hydrogen. 

The asymptotic limit of the radial continuum wave- 
function of the final state with angular momentum L is of 
the generic form [24| : 

Nk 
r 



lim Rkxir) 



- s\n[kr + $fe^L(f)] 



(11) 



where N}; — ^2/(7rfc) is the normalization constant in the 
energy scale and the phase has the general form: 



$fe,L(r) = Z\n{2kr)/k + T^L{k) - 



(12) 



We note that this phase includes the logarithmic diver- 
gence characteristic of the Coulomb potential of the ionic 
core with charge Z, in the asymptotic region. The Coulomb 
potential influences also the scattering phase-shift rjL{k), 
which can be rewritten under the form: r]L(k) — ^^{k) -\- 
5]^{k) where (7^ = arg[r(L -f 1 — iZ/k)] is the Coulomb 
phase-shift and where the correction S^i^k) originates from 
the short range deviation of the ionic potential from a pure 
Coulomb potential, see for instance [211 [21]. Obviously, in 
the case of an hydrogenic system, one has ^^(fc) = 0. 

To derive the asymptotic form of the perturbed wave- 
function Pk,a('")j it is in principle enough to establish the 
limiting structure of the differential equation it verifies. 
From Eq. one observes that, in the asymptotic limit 
r — >■ cx), the second member vanishes, as a result of the ex- 
ponential decay of the bound state wavefunction R-mJi ir). 
One is left with a standard Schrodinger equation for posi- 
tive energy which is solved by imposing outgoing wave 
boundary conditions to the solutions ^32j : 

lim pK,A(r) oc — exp [i{Kr + $«,>(?■))] . (13) 

r—^oo r 



It is also of interest to derive explicitly the limiting forms 
of the terms entering the expression of pK,x{'r) in Eq. ( 10 ). 



Regarding its real part, the sum over the discrete states v 
can be neglected, as each term goes asymptotically to zero. 
Thus, for large r, it reduces to a principal part integral: 

i.ik.Wl^W'^.- ^--'^"^--"'!^-''' . (14) 

This integral can be estimated by extending the integra- 
tion range e^' —oo and replacing the radial continuum 
function Rk',x{'>') by its asymptotic limit according to the 
prescription in Eq. (11). Then, writing the sine function 



under its exponential form and performing contour inte- 
grations, with semi-circles around the pole at e^, one gets: 



lim ^W,x{r)] 



-cos[Kr + ^,^x{r)]{Rn,x\r\R„ 



(15) 



Regarding the imaginary part given by the last term in 



Eq. (10), it is enough to substitute again the asymptotic 



form of RK.xir), so that: 
lim 5[pK.A(r)] w -^^^sm[Kr + ^K.x{r)]{RK,x\r\RniA)- 

r— foo r 

(16) 

Then by regrouping the real and imaginary parts, one gets 
the final expression: 



lim [PkA'^) 



ttN, 



- exp [inr + i^n,xir)] (i?K,Ak|i?„,,£.), 

(17) 



which corresponds to a complex outgoing wave [32], as 
expected from Eq. (13), weighted by the dipole matrix 
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element associated to the one-photon transition from the 
initial state. We note that adopting the so-called "pole- 
approximation", which consists in neglecting the off-shell 



part {i.e. the real part given in Eq. (15)), would lead to 
a loss of the phase information of the process since the 
perturbed wavefunction then would be a standing wave 
rather than an outgoing wave. 

The corresponding asymptotic approximation for the 
second-order radial matrix element, Eq. ([6]), is obtained by 
substituting in Eq. ([s]) the asymptotic expressions, Eqs. ( 11 ) 
and ( 17 1, for the radial wavefunctions of the final state and 



of the intermediate state, respectively. One has explicitly: 

/>oo 

X / drsin[/cr + $fe,z,(r)] r e['('"'^+*'''^(''»l. 
Jo 

(18) 

To introduce the next step in our approximate treatment, 
one rewrites the sine in its exponential form and develop 
the expressions of the phases ^K..\{r) as given in Eq. (12). 
One is left with two distinct contributions containing in- 
tegrals either of the type J+ or J_ that are defined as 
follows: 

J± = ±-j dr r^+^z{i/K±i/k) j.^^ ^ f^y^ 



± 



1 



2i \ K ± k 



2+iZ{l/K±l/k) 



T[2 + iZ{l/K±l/k)], 

(19) 



where we have used an integral representation of a Gamma 
function T{z) with complex argument. In our case, the 
contribution of the J+ integral is vanishingly small as com- 
pared to that of J_. This is due to the IR photon energy 
being small compared to the final kinetic energy of the 
electron, cj — fc^/2 — k^/2 <^ fc^/2, so that the difference 
|k — fc| w uj/k is much smaller than the sum k + fc « 2fc ± 
oj/k. As a result, the fast oscillations of exp[i{K + k)r] lead 
to a relative cancellation of the corresponding integral, as 
compared to the one containing the factor exp[i(K — k)r]. 
Neglecting the J+ contribution, the asymptotic expression 
reduces to: 

TLaA(fc;e«) « J NkN^ (i?«,Ak|i?„.A) 
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exp 
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4)] 



|k-/c|2 

(2k)*^/'= r[2 + iZ{l/Hi - 1/fc)] 
^ (2fc)*^/fc {k - kyz{i/K-i/k} 



(20) 



which was used by us |3j to obtain estimates of the phases 
occurring in two-photon transitions entering RABBIT tran- 
sition amplitudes. The first two lines in Eq. (20) are real. 



they contain a one-photon matrix element from the bound 
state into the continuum, but also an exponential factor 



describing the strength of the continuum-continuum tran- 
sition from K to fc. The exponential factor decreases with 
the probe photon energy, uj = k'^ /2 — /2, which indicates 
that large energy leaps in the continuum are strongly sup- 
pressed. At a given laser probe energy, however, the expo- 
nential factor increases with the final momentum, fc, which 
indicates that it becomes easier for the photoelectron to 
interact with the probe field. The third line is a simple 
phase factor containing the scattering phases of the con- 
tinuum states. Finally, the fourth line is a complex factor 
that depends on three quantities: final momentum, fc; the 
laser probe frequency, uj; and the charge of the ion, Z. 

A more formal derivation of this result, based on a 
closed-form representation of the Coulomb Green's func- 
tion is given in the Appendix A. [We have found a typo in 
Eq. (7) in ref. [3^: The ratio, (2fc)*/V(2fca)*/'=° , should be 



inverted, as is evident from Eq. (20) in the present work]. 

Thus, in the asymptotic limit, the phase of the radial 
component takes the form: 

arg [rL,AA(fc;e«)]«|(i-A-i) 

+ V\(K)-VL{k) + (f>cc{k,K), (21) 



where 



0cc(fc, k) 



arg 



(2k)*^/« T[2 + iZ{l/K-l/k)] 
(2fcpA {k- kyz{i/K-i/k) 



, (22) 



is the phase associated to a continuum-continuum radia- 
tive transition resulting from the absorption of w, in the 
presence of the Coulomb potential, Z. It is important to 
note that it is independent from the characteristics of the 
initial atomic state, as well as from the amplitude of the 
field. It is illustrative to study the continuum-continuum 
phase in the limit of a small photon energy, w « fc(fc — k), 
which yields a simplified expression: 



0(r^*)(fc;c.) 



arg 



2fc2 

UJ 



iZuj/k^ 



T[2 + iZuj/k^ 



(23) 



where it becomes clear that it is the product: Zui/k'^, 
which determines the size of the continuum-continuum 
phase. This expression is expected to be valid in the so- 
called "soft-photon" limit, k'^/2 ^ uj, where the exchange 
of energy w and the corresponding momentum transfer 
Afc = uj/c do not significantly modify the electron state 
[23]. The substitution, uj —uj, yields the (Fourier) 
component corresponding to stimulated emission of light. 
We note that the phases corresponding to absorption and 
emission have opposite signs, but that they are otherwise 
identical in the soft-photon limit. 



Replacing the formula obtained in Eq. ( 20 ) for the ra- 
dial component in the expression of the full transition am- 
phtude M{k,ei + il) given in Eq. (jsl, one gets its general 



6 



form in the asymptotic limit 
2^2 



1 



exp 



ttZ / I 1 

~2~ \ k~ k 



{2Ky^/^ r[2 + iz{i/K - i/k)] 

^ (2fc)'^/'= (k - /c)»2(i/«-i/'=) 

(24) 

To address the question of its phase, one notices that be- 
sides a trivial contribution from the spherical harmonic in 
the final state, YL,mi(fc)j it contains only phase-shifts that 
are governed by the angular momentum A of the intermedi- 
ate state, i.e. a state that can be reached via single-photon 
ionization. More precisely, for a given transition channel 
characterized by the angular momenta of the intermedi- 
ate and final state £i —?' X ^ L , the phase of the matrix 
element reduces to: 



arg [ML,xji (fc, e^)] = TT + arg[YL,„, (fc)] + (/)o + 0^ 



(25) 



where (pQ and (pi^ are the phases of the XUV field fl and 
of the IR laser w, respectively. We stress that the final 
state scattering phase, 7]L{k), cancels out and that it enters 



neither in Eq. (24) nor in Eq. (25). 



Eq. (25 ) represents one of the major results of our the- 
oretical analysis. It shows that, within the asymptotic 
approximation and besides trivial spherical harmonic con- 



tributions and the phases of the fields [line 1 in Eq. ( 25 



the phase of a two-color ATI transition amplitude has two 
components: i) One is directly linked to the quantum- 
mechanical phase-shift of the one-photon XUV ionization 
amplitude, here — 7rA/2 -t- rixin); ii) The other, denoted 
(pccik^n), is in some sense "universal", it describes the 
phase brought by the absorption of the probe photon a;, 
in the presence of the Coulomb potential with charge Z . 
Then, as shown below, when comparing laser-assisted ion- 
ization originating from distinct atomic states, the energy 
derivative of the phase-shifts in the first term contribute 
to a Wigner-like time-delay. On the other hand, the dif- 
ference between the "universal" terms 0cc gives rise to 
a measurement-induced delay, associated to continuum- 
continuum stimulated radiative transitions in the presence 
of the Coulomb potential of the ionic core. 

In Fig.[3j we present the continuum-continuum phases, 
associated with absorption (red) and emission (blue) of 
a probe photon leading to the same final energy, calcu- 
lated using the asymptotic approximation, Eq. ( |22[ ). The 
continuum-continuum phases for probe photon absorption 
are positive while those for stimulated emission are neg- 
ative, but approximately equal in absolute value. In the 
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Figure 3: Continuum— continuum phases for absorption (red, up- 
per curves) and emission (blue, lower curves) calculated using the 
asymptotic approximation (dashed curves) [Eq. |22] and the long- 
range amplitude-corrected asymptotic approximation (full curves) 
[Eq. |30| . These approximate phases are compared with the exact 
calculations (black symbols) from the Is state in hydrogen. The ex- 
act results for final angular momentum L = {+ symbol) and L = 2 
(x symbol) are computed by subtracting the one-photon scattering 
phase from that of the exact two-photon matrix elements. The data 
correspond to Z = 1 and to a laser probe with uj = 1.55 eV. It is 
interpolated between the discrete harmonic orders. 



next subsection, we show that, still in a single- active elec- 
tron picture, it is feasible to improve the accuracy of our 
approximate treatment with the help of semi-classical ar- 
guments. 

3.3. Long-range amplitude effects 

In order to go to the next level of our asymptotic ap- 
proximation, we must include not only long-range phase 
variations of the continuum states, but also long-range 
variations of the amplitudes. Indeed, the normalization 
constants contained in the asymptotic forms of the radial 
functions Rk,L{i') and Pn,\{r) can be modified to account 
for the long-range influence of the Coulomb potential. For 
instance, the modified final state normalization constant 
is: 



Nkir) 



iTp{r) ' 



where 



p{r) = ^2{e-V{r)) 



V{r)/k, 



(26) 



(27) 



is the local momentum from Wentzel-Kramers-Brillouin 
(WKB) theory [33]. The same remark applies to N^{r) 
for the perturbed radial function. Long-range amplitude 
effects can be approximated by expanding the quantities 
Nk{r) and N^fr), 



Nk{r)N,{r) 



4 

n'^kn 



1 



1 



1 

fc2 



(28) 
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to first-order in the Coulomb potential. The second term 



within the brackets in Eq. ( 28 ) contains the first-order am- 



plitude correction to the matrix element. Evaluation of the 
long-range amplitude contribution leads to a correction to 
the continuum-continuum phase: 



arg 



1 



iZ 
T 



1 



1 

fc2 



l+iZ{l/K- 1/fc) 



The final continuum-continuum transition phase is: 



(29) 



(30) 



where <j)cc{k, k) is given in Eq. ([22]). In Fig. [3] we show that 
including such long-range effects improves the accuracy of 
the approximation greatly, leading to accurate quantita- 
tive result already at relatively low energies in the contin- 
uum. The accuracy in the lower energy range can be fur- 
ther improved by using a regularization method designed 



to remove the effect of the singularity in Eq. (28) as r, k 
and K go to zero. 

Going beyond the approximations given here, namely 
performing exact ah initio computations of the matrix el- 
ements M(k, ti ~\- 17) in polyelectronic systems, is out of 
reach of present computational capabilities. It is only 
in the special case of hydrogenic systems, that such 2- 
photon amplitudes can be computed with arbitrary preci- 
sion. Thus, with objective to delineate the range of validity 
of our approximation, we turn now to a brief presentation 
of the "exact" calculations in hydrogen. 

3. 4- Exact calculations of 2-photon ATI matrix elements 
in hydrogenic systems 
The principle of the calculation is outlined here for s- 
states. Numerical data for other states will be given below. 
We first express the transition amplitude given in Eq. ^ 
for the case £i — 0, m; = which implies A = 1 so that the 
angular momentum of the photoelectron is either L = 0,2. 
Accordingly, two distinct amplitudes contribute to ATI 
transitions like the one depicted in Fig. [l] (c): 



M{k,e,: 



3i 



4e'"'^'^^2,o(^)T2,i,o(fc;e,) 



(31) 



where the radial components for s— states are of the form: 

TL,i,o(fc; ^^)k^o = {Rk,L\r\Gi{e^)\r\Rn,,o}, (32) 

with L — 0, 2; Ui labels the initial atomic s— state and 
G\=i{£k) is the radial component of the Coulomb Green's 
function for angular momentum A — 1. The general form 
of the Green's function with energy argument is: 

GA(r ,r;e«, = hm V — — . 33) 

e^O+jH - El/ + te 



As already mentioned, the infinite sum over the index v 
runs over the whole (discrete + continuous) spectrum of 
the hydrogenic system. Closed form expressions for Gx 
are known, see for instance [33| . Here, we have used the 
expression given as an expansion over a discrete Sturmian 
basis: 



Gx{r\r;ei^) = 



E 



Si^^X,x {r')S^^x,x{r) 



1 



(34) 



where x — \/ —^e^, and the so-called Sturmian functions 
Sv,x,xij'^ have a structure similar to the bound-state hy- 
drogenic radial functions [331 '■ 



S„.x.x{r) = 2xa 



X e--^2xr)^Llt\l,i2xr), 



(35) 



where L1^^_^{z) are associated Laguerre polynomials. In 
the amplitudes for ATI transitions, = ei + ft > 0, 
and the quantity, x — i^/2\e^\, is pure imaginary. It is 
then convenient to use Pade-like resummation techniques 
to compute the infinite sum over the index see, for in- 
stance ref. [31]. 

In Fig. [3] we present the exact continuum-continuum 
phases from the Is state in hydrogen. These phases are 
defined as the total phase of the exact matrix element, 
-^L,i,o(fc; £k)j minus the one-photon phase [see line 2 of 
Eq. (25)]. Our approximate calculation including long- 



range amplitude effects, Eq. (30 1, is in excellent agreement 



with the exact calculations except at low energy. 



3.5. Phase of the classical dipole 

Finally, we present a simplified derivation of the continuum- 
continuum phase, 4>cc(k,K), using a classical approach. 
The dipole associated with the absorption of radiation at 
frequency w by a free electron in the presence of a Coulomb 
potential, can be calculated using Larmor's formula, 

/•OO 

dc{k;uj) — / dt rk{t) cxp[—iujt], (36) 
Jo 

where it is assumed that the electron follows a field-free 
trajectory, rk{t), that starts close to the ion, r'fc(O) « 0, 
and then moves out away from the ion with an asymptotic 
velocity, k. The integral can be cast from time to space 
using the r— dependence of the velocity: 



Vk{r) = _ 2V{r), 



(37) 



where fc^/2 is the final kinetic energy of the electron at 
large distance from the ion. Using the differential dt = 
dr/v{r), the time can be written as i(fc; r) = dr' /vk{r') + 
C, where C is an integration constant. In the case of the 
Coulomb potential, V{r) = —Z/r, the integral becomes 



i(fc;r) 



dr' 



^fc2 + 2Z/r' 



+ C 



(38) 
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in the asymptotic limit, i.e. when ^ Z/r. This 

provides an approximate time-position relation valid at 
large distances from the origin. In the special case where 
the electron starts from the origin [t, r] = [0,0], the ex- 
act integration in Eq. Q leads to C = - Z \a[2k'^ / Z] / . 
Keeping for the moment this value of the constant and re- 
placing the asymptotic form of the time in the expression 



of the dipole Eq. ( 36 ) , one gets 



dc{k]uj) ~ J dr rexp 
k Jo 



1 /2k^ 



2 \ iZuj/k „co 



k \ Z 
k 



dr r^+^z^/'^' 



exp 



exp 



/o 
3t:Zuj 
2k^ 



-i—r 
k . 



2fc^ 
Zu 



3 s iZu/k" 



T{2 + iZu/k^), 



(39) 



where the next-to-last line is real, with an exponential fac- 
tor that decreases with w, but increases with k, in excellent 



agreement with the quantum mechanical result, Eq. (24) 



Furthermore, the last line contains the complex gamma 
function times an algebraic factor, also in close connection 
to the quantum counterpart. Clearly, the dipole corre- 
sponding to absorption is a complex quantity with phase: 



0c(fc; 



arg [dcik;uj)] 
arg 



. (40) 



which is closely related, but not identical to its soft-photon 
quantum counter part , given in Eq. ( 23 ) . 



In the quantum mechanical case, the electron starts in 
a bound state with some spatial extent and not exactly 
from r = 0. In order to account for this uncertainty on 
the initial position, we may choose a different value of the 



integration constant C, introduced in Eq. (381, in order 
to come closer to the quantum mechanical dipole. This 
matching-procedure is reminiscent of the method used in 
ref. [15], to determine the "best" initial radial position 
for the electron within the eikonal Volkov approximation. 



Within the lowest-order approximation of Eq. ( 38 ), we find 



a simple relation between the initial position and the in- 
tegration constant: 



ro « exp 



Ck^ 



(41) 



which is valid at high-energy. It is convenient to set C = 
— Zln[2k]/k^, a choice corresponding to an initial radial 
position To ~ 1/(2A:), as was identified in ref. [12 . Clearly, 



phase if we want to deduce the "true" initial position. In- 
terestingly, in our approach this exact inital position is 
not critical. In fact, our results are stable with respect to 
rather substantial modifications of the wavefunctions close 
to the core. It is the behavior of the wavefunctions far away 
from the core that must be described accurately, using the 
asymptotic approximation. We now turn our attention to 
the applications of the complex ATI matrix elements, to 
the determination of attosecond delays in photoionization. 

4. Attosecond Time-delays 

The complex amplitude M (k, + O) (denoted M^"^ 
for conciseness in the following) for the joint absorption 
of an XUV photon and of an IR laser photon w, is of 
course not a direct observable in any experiment. Only 
the square modulus of a complex transition amplitude can 
be measured. Thus, if M^°-'> is the only amplitude leading 
to a given final state, there is no way to determine its 
absolute phase. However, in the cases of interest here, 
with attosecond XUV pulses in the presence of a probe IR 
laser field, several other channels are open which can lead 
to the same final state, thus making it feasible to observe 
phase-dependent interference patterns. This property is 
exploited in the following schemes. 

4.1. Attosecond delay measurements using pulse trains 

In the RABBIT scheme there are two dominant com- 
plex amplitudes (quantum paths), Af(°) and Af^*^', asso- 
ciated with the absorption of harmonic H2q-i or H2q+i 
plus absorption or emission of a laser photon with phase, 
ibwT = ±(/);^, leading to the same final sideband, S2q- The 
photoelectron will transit via different intermediate states: 
K< and K>, corresponding to different intermediate ener- 
gies: e< = e,; + 2quj ~ ijj — k\/2 and e> = -I- 2gcj + uj = 
K^/2. Clearly, the energy of the corresponding intermedi- 
ate states are one photon below and above the sideband, 
as is illustrated in Fig. [l] (a) . This implies that the mea- 
sured intensity of the sideband, S2q, will depend on the 
phase difference between the two quantum paths: 



when the first-order amplitude correction is included, Eq. (30 1 
the initial position is adjusted accordingly. The continuum- 
continuum phase being only one part of the total quantum 
mechanical phase, we have also to include the scattering 



+ 2|A//('')||Af('=)| cos arg(A//(")*Af('=)) (42) 

which corresponds to a standard interference phenomenon, 
as summarized in Eqs. ([T]) or, equivalently in ([2| displayed 
in the introduction (Sec. [T]). If we now apply Eq. (25), 
and assume that the total contribution to the complex 
amplitudes can be approximated by a single intermediate 
angular momentum component. A, the phase difference is 



arg 



M 



+ V\ii^>) - V\{i^<) + (t)cc{k,Ky) - cj)cc{k,K<^) . (43) 
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Here the first line contains the phases of the fields, while 
the second line is in terms of scattering and continuum- 

continuum phases. The corresponding measurement delay, 
Ts as defined in Eq. ([2]), is the sum of finite-difference ap- 
proximations to a Wigner-like time-delay: 



??a(k>) ~ rix{i^<) 
2w 



and to a continuum-continuum delay: 



„(fc;a;) 



2w 



(44) 



(45) 



so that Tg T\ + Tec- 

In a more general case, the interference will result from 
complex amplitudes with two terms for absorption and for 
emission, M^"/'^) = M(''/'''+) -I- M^"/'' corresponding to 
different intermediate angular momenta, X = £i ± 1. The 
interesting point is that the Tcc will not be affected, since 
the phases (f>cc are independent of A as well as of the fi- 
nal angular momenta L = £i,£i ±2, see Eq. (22). The 



remaining part can be interpreted as an effective Wigner- 
like delay, f\^i.±i, which has to be computed taking into 
account the relative weights of the four components enter- 
ing the expression of the transition probability amplitude. 

4-2. Attosecond delay measurements using single pulses 

In this section we apply the perturbative treatement to 
laser-assisted photoionization by a single-attosecond pulse 
(SAP) and a probing laser field. The ionizing attosecond 
field is 



E. 



SAP 



it) 



dn Ene^p[~int]/2'!T, 



(46) 



with Efi — \Eq\ exTp[i(j)fi\ being complex Fourier coeffi- 
cients. While our approach is very general, it can be il- 
lustrated by considering a Gaussian frequency distribu- 
tion centered on the frequency fio, as depicted in the left 
panel of Fig. |4] The probing laser field is assumed to be 
monochromatic and real, E{t) = 2\E^\cos[ijj{t — r)], so 
that it can account for both absorption and emission pro- 
cesses. Our perturbative approach cannot be used to fully 
account for the large momentum shifts that are typical of 
streaking spectrograms, but it does allow us to study the 
onset of streaking. As we shall see below, streaking-like 
behavior is clear already in the perturbative regime. In the 
right panel of Fig. |4j we sketch the interfering processes 
to a certain energy e: (d), (a) and (e) from lowest-order 
perturbation theory. The dominant contribution (d) cor- 
responds to absorption of a single XUV photon, while the 
upshifted (a) and downshifted (e) photoelectron spectra 
corresponds to absorption and emission of an additional 
laser photon cu, respectively. 

The matrix element for one-photon ionization with a 
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Figure 4: Sketch of the quantum paths describing the onset of 
"streaking" for a photoelectron ionized by a single attosecond pulse 
and probed by a monochromatic laser field. The first-order pho- 
toelectron wave packet is centered at to = ti + f!o- Within the 
bandwidth of the attosecond pulse, any energy, e, can be reached 
by path (d), where a single XUV photon with frequency f! is ab- 
sorbed. Alternatively, the same energy can reached by path (a) by 
absorbing a less energetic photon, f2<, and a laser photon, oj; or 
by path (e) by absorbing a more energetic photon, f!>, and then 
emitting a laser photon, —u). At the high-energy end of the photo- 
electron distribution (indicated by e) , the dominant contributions to 
the photoelectron wave packet are (d) and (a), due to f!> being far 
in the upper energy range of the XUV bandwidth. 



photoelectron emitted along the polarization axis is 
^3/2 ^ . /2A + 1 



M('^)(H;e, + 17) = ^^£;oV(- 

7 Z ^ 



irixik) 



X {Yx,o\Yi.o\Yi,,o){Rk.x\r\Rn,A) ^m^.O, 

(47) 

where the XUV frequency satisfies ei + fl = e = 12 
for process (d). Note that only initial states with zero 
magnetic quantum number, rrii — 0, will contribute to 
photoelectron emission along i, as indicated by the Kro- 



necker delta at the end of Eq. (47) that originates from 



the explicit properties of the spherical harmonics along z: 
Fa,™. (5) = x/(2A + l)/3 

The two-photon matrix element for photoelectrons along 
the polarization direction is 



(Stt) 



3/2 



2L + 1 



M{kz-t, + n^)='-^^^^En^E^Y.^-') V -x 

% L ^ 

X Y.(Yl.o\Yi,o\Yx,o){Yx,o\Yi,o\Yu,o) 

X 

X TL^A,f,(fc; Q + ri<) 6,ni,0, 



ir]L (fe) 



(48) 



where ei-\-rt^-\-uj = e for process (a). A similar expression, 
M{kz, ti -f r2>), can be written for process (e) where -|- 
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S7> — cj = e, and where the complex conjugate of the laser 
field, E^, is used for emission of a laser photon. Again, 
only initial states with zero magnetic quantum number 
contribute. 

The PROOF method [20] relies on the interference be- 
tween these three quantum paths, denoted M'^'^\ M^"' and 
M'''^\ as a function of the XUV-laser delay, r. The prob- 
ability for emission of an electron with final momentum k 
along the polarization axis z is 



P{kz) (X |M('*) + M'^"^ + M(")p 



(49) 



where only cross-terms lead to delay dependent modula- 
tions. Provided that the laser field is relatively weak, we 
expect the interference cross-terms (d)-(a) and (d)-(e), will 
dominate over the cross-term (a)-(e), because the former 
involve the exchange of only one laser photon while the 
latter contributes to higher order terms in a perturbative 
treatment, involving the exchange of two laser photons w. 

In order to illustrate the origin of the attosecond delay 
within this framework, we consider a final energy that is 
higher than the central energy of the photoelectron wave 
packet, e = ti + ^l > ti + ^Q. In this region, the uj modula- 
tion is dominated by the (d)-(a) interference, because path 
(a) goes through a central part of the XUV distribution, 
while path (e) passes through a higher frequency range, in 
the upper part of the XUV bandwidth [Fig.|4]. Increasing 
the bandwidth of the attosecond pulse or decreasing the 
probe photon frequency makes this distinction less pro- 
nounced. We restrict our analysis to the case of an initial 
s-state, li = 0, and study the interference of the cross term 
(d)-(a). The relevant phase reads 



arg 



(50) 

where we have used the asymptotic approximation for the 
phase of the two-photon matrix element as in Eq. ( 25 ) , 



thereby, introducing the continuum-continuum phase into 
the framework of laser-assisted photoioization by single 
attosecond pulses. Also, we note that the phase-shifts 
present in the two-photon matrix element are those of 
p— waves (A = 1). In terms of the temporal delays, the 
UJ modulation is displaced by 



m + Tx{k) 



UJ 

Tccik;uj) 



TT 

2^' 



2uj 
(51) 



the electron distribution using the cross term (d)-(e), and 
we find the same result but shifted by a half laser period, 
i.e. out of phase by tt. This tt— shift of the modulation 
is important to explain the onset of streaking, because it 
ensures that the high-energy probability of the electron 
spectra is maximized when the low-energy part is mini- 
mized. In the central region, e « -I- VIq, the (d)-(a) and 
(d)-(e) contributions will be comparable, leading to a rela- 
tive cancellation of the uj modulation, leaving only the 2uj 
modulation. In total, the electron momentum distribution 
is slightly shifted up or down depending on the sub-cycle 
delay between the laser field and the attosecond pulse: It 
is streaked as expected from the classical (or strong-field) 
picture. Using Eq. (51), we find that this streaking mod- 



ulation is displaced by a Wigner-hke delay, T\{k), and by 
the continuum-continuum delay, Tcc{k]uj). We can inter- 
pret this delay as the time it takes for the electron to be 
photoionized plus the time it takes for the measurement 
process to occur, i.e. for a lower energy electron to ab- 
sorb one probe photon so that it may interfere with the 
direct path. We stress that this analysis was made as- 
suming a monochromatic probe field, while it is common 
in streaking experiments to use few-cycle IR laser pulse, 
which leads to a convolution (a blurring effect) of the mo- 
mentum modulation and of the attosecond time-delays. 

4-3. Comparison between the two measurements 

We now briefly discuss the small differences of Tcc oc- 
curing between the RABBIT method and the streaking 



method. Eq. (52) implies that 



Tcc{k\Uj) 



0cc\ f^j ^> I H^cc 



2uj 



(f>ccik,K>) 



OJ 



(53) 



where we have used the definition of the continuum-continuum 



delays Tcc, Eq. (45), and the following relation for the 



continuum-continuum phases 



(f>ccik,K<) W -(t>cc{k,Ky), 



(52) 



which is exact in the soft-photon limit. Similarly to Eq. ( 51 ) 



we can compute the modulation at the low energy end of 



where the first line corresponds to the RABBIT method 
and the second line corresponds to the streaking method. 
These different "flavours" of Tcc merge completely in the 
soft-photon limit, e ^ uj, but differ slightly at low kinetic 
energies, where the phases are not exact opposites. In 
both methods, Tcc varies with k inside the bandwidth of 
the photoelectron wavepacket. 

Furthermore, it is interesting to note that the group de- 
lays of the attosecond pulses and Wigner delays of the pho- 
toelectrons appear as discrete derivatives over one and two 
OJ photons in the streaking and RABITT methods respec- 
tively. Here, we have assumed that all delays are slowly 
varying so that the discrete derivatives are equivalent. We 
also note that the streaking method relies on the cross- 
terms: (d)-(a) and (d)-(e); while the RABBIT method re- 
lies on the (a)-(e) contribution. The weaker signal in the 
RABBIT method is not a problem, since it is recorded be- 
tween the harmonics, on "zero background" . Consequenty, 
the modulation frequency of the streaking/PROOF signal 
is UJ due to the single laser photon involved in each ap- 
propriate cross term, while the corresponding RABBIT 
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frequency is 2a; due to the two laser photons involved in 
the latter cross term. Another important difference be- 
tween the methods is related to symmetry (parity): The 
w signal is somewhat restricted to the z direction and it 
is not observed if photoelectrons are collected in all direc- 
tions; while the 2a; signal is more general and it prevails 
also when electrons in all directions are collected. Inter- 
estingly, we have shown that all three attosecond char- 
acterization methods (RABITT, PROOF and streaking) 
provide equivalent temporal information about the XUV 
ionization process, even though they are built from differ- 
ent sets of cross-terms associated to different interfering 
quantum paths. 

5. Results 

5.1. Calculations of attosecond time-delays 

In Fig. [5j we plot atomic delays relevant for laser- 
assisted photoionization. The Wigner-like delays, t\, are 
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Figure 5; Atomic delays of laser-assisted photoionization in hydro- 
gen for three different angular momenta of the wave-packet wave- 
function in the intermediate state (s, p and d : A = 0, 1 and 2). The 
delays are calculated using the regularized approximation. The data 
correspond to Z = I and a laser probe with u) = 1.55 eV. 

calculated from the finite-difference derivative of the scat- 
tering phase of hydrogen for angular momentum, A = 
[0, 1, 2], corresponding to s, p and d continuum wave pack- 
ets. The universal continuum-continuum delay, Tcc, plus 
these Wigner-like delays yields the total atomic delay, Tg. 
Notice that the sign of t\ and t^c are opposite, so that 
Te is smaller than either of the contributions individually. 
The increase of t\ with the angular momentum can be 
understood as due to the repulsive, short-range, centrifu- 
gal potential. The total delay is negative in this example, 
which implies that the electron appears as being advanced 
compared to the probe field. 

In Fig.[6l we evaluate the accuracy of Tcc as derived by 
various degrees of the asymptotic approximation by com- 
parison with exact calculations in hydrogen from the Is 



state along two different angular momentum sequences: 
s — > p — > s (+) and s — > p — > d (x). Indeed, we find 
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Figure 



Continuum-continuum delays calculated using the 



asymptotic approximation (dashed curve), Eq. l |22| , the long- 
range amplitude-corrected asymptotic approximation (thick curve), 
Eq. l |30| , and the regularized asymptotic approximation (dash-dot 
curve) . The approximate delays are compared with exact delays 
computed from the Is state in hydrogen by subtracting the interme- 
diate Wigner delay. The data correspond to Z = 1 and to a laser 
probe with u = 1.55 eV. 

that the exact Tcc are almost completely independent of 
the final state angular momentum. The asymptotic ap- 
proximation (dashed curve) predicts the correct qualita- 
tive behavior of the delay, but it slightly overestimates its 
magnitude. By taking into account the long-range am- 
plitude effects, the agreement is excellent at high kinetic 
energies. The disagreement at low energy can be removed 



by avoiding the radial singularity in Eq. (28), namely by 
an ad hoc substitution of r — >■ r -I- j(l — \k ~ k\/2). In this 
way, we obtain a "regularized" continuum-continuum de- 
lay (dot-dashed curve), which is excellent at all energies 
in the range. 

In Fig. [7) we address the question of the universality of 
Tcc by computing the exact two-photon phases from three 
different intial states in hydrogen: Is, 2s and 2p. We 
find that all seven different angular-momentum sequences, 
£i X ^ L, line up on the same universal curve, in excel- 
lent agreement with the regularized continuum-continuum 
delay. In this way, we have verified that not only s— type 
initial states have similar Tcc, but also initial p— states with 
non-zero angular momentum. 

In Fig. [Sj we present a contour plot that provides a 
rough overview of the magnitude of Tcc as a function of 
the kinetic energy of the photoelectron and of the wave- 
length of the probe field. In general, Tcc decreases with the 
kinetic energy and it increases with the wavelength of the 
probe field. We note that a softer probe photon leads to an 
increased delay. A few selected Tcc curves are displayed in 
Fig. |9] for a more quantitative comparison at some exper- 
imentally relevant wavelengths. We conclude that Tcc is 
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Figure 7: Continuum-continuum delays from exact calculations in 
hydrogen from different intial states: Is (black symbols), 2s (red 
symbols) and 2p (blue symbols). This data demonstrates the uni- 
versality of the continuum— continuum delay and that it is valid not 
only for initial states of s— character. The exact data is in excellent 
agreement with the regularized approximation. The data correspond 
to Z = 1 and a laser probe with uj = 1.55 eV. 



not extremely sensitive to the probe wavelength, except at 
low kinetic energies, where it takes larger negative values 
when the probe wavelength increases. As already men- 
tioned, the delays converge to zero as the kinetic energy 
is increased, however, the convergence is rather slow and 
there is still ~ —10 as of delay remaining at a relatively 
high kinetic energy of ^ 100 eV. 



Finally, in Fig. 10 we examine the validity of the soft- 



photon approximation by comparing r^c and ric"'^*'* , calcu- 



lated using Eq. ( 22 ) and Eq. ( 23 ) , respectively. The soft- 



photon limit over-estimates the magnitude of the delay at 
low energies, but it converges rapidly towards the quantum 
mechanical result as the kinetic energy is increased. 

5.2. Conclusions 

The main result of this work is the determination of 
the phase associated with two-photon transition matrix 
elements for laser-assisted XUV photoionization, namely 
Eq. (25). These phases have broad apphcations in at- 



tosecond science as they enter naturally in most charac- 
terization methods, as well as in quantum control schemes 
and delay experiments of more general character. We 
have demonstrated that the phase is composed of two dis- 
tinct atomic contributions: i) The one-photon scattering 
phase of the intermediate state and ii) a quantity, given 



in Eq. (22 1, that we call the continuum-continuum phase. 



The latter is universal and it describes the added phase 
induced by the transition from the intermediate contin- 
uum state to the final continuum state. It is independent 
of the short-range behavior of the atomic potential and it 
depends only on three quantities: the final momentum, the 
laser-probe frequency and the charge of the remaining ionic 
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Figure 8: Contour plot of continuum— continuum delays at different 
kinetic energy and laser probe wavelength. The delays are calculated 
using the regularized asymptotic approximation. 



core. As expected from the strong-field approach, the tran- 
sition phase is reduced to that of the one-photon (interme- 
diate) scattering phase when the charge of the remaining 
ion is neglected, i.e. for a short-range potential. Other 
interesting findings are related to the classical-quantum 
correspondence: First, the phase of the classical dipole, 
for exchange of radiation in the continuum, Eq. (40), is 



closely related, but not identical, to the soft-photon limit 
of the quantum mechanical continuum-continuum phase, 
Eq. (23). Further, as epitomized in Eq. (41), there ex- 



ists a relationship between the initial radial position of 
the ejected electron (a classical concept) and the asymp- 
totic quantum phase. We have also demonstrated how the 
theory of complex transition matrix elements, originally 
developed for RABBIT, can help to better understand at- 
tosecond streaking measurements. In our interferometric 
interpretation of streaking, we stressed the benefits of a 
weak and monochromatic laser probe field. Given such ex- 
perimental conditions, quantitative analysis of more com- 
plex systems can be carried out using two-photon transi- 
tion matrix elements, corrected by many-body perturba- 
tion theory. 
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Appendix A: Formal derivation of the approximate 
expression for the radial component of the 2-photon 
ATI transition amplitudes 

The radial component of the two-photon ATI matrix 
elements given in Eq. can be rewritten: 

TL,x,i,{k;e, + n) = lim {Rk.L\r Gx{ei + n + ie) r\Rr,^^t^^). 

(54) 

Here the initial atomic state is represented by the ra- 
dial wavefunction i?„^ £. (r), while Rk,L{f) is the wavefunc- 
tion of the photoelectron with energy k"^ /2 = Ci + fl ± uj. 
Ga(?'2,?'i; Ei + ^) is the radial component, of the Green's 



function given in Eq. ( 33 1 for angular momentum A and 
positive energy argument = -|- O > 0. We note the 
presence of the positive imaginary infinitesimal ie in the 
argument of G: It corresponds to the change k — ^ n + ie/K, 
which ensures the presence of a converging factor e^'^^ in 
the integrals below, which otherwise would be divergent. 
For the sake of conciseness, we will omit this factor in the 
following. 

In the limit of large values of the coordinates ri and 
r2, the Green's function takes the following limiting form, 
written for the case of a Coulomb potential, [371 EH] : 



G\{r2,ri;ei + fl) ^ ' — e 

rir2K 



i[Kr>+$^_;^(r>)] 

X sin[Kr< + $K,A (?•<)] (55) 



where r> (resp. r<) denotes the larger (resp. smaller) of 
(ri,r2) and the phases <&«.>(?') are defined in Eq. (12). 



We consider now the case of two-photon ATI from the 
Is state of a hydrogenic system with Ri o{r) — Ci.oe"^'" 
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Figure 10: Continuum-continuum delays calculated using the 
asymptotic approximation (dashed curve) using Eq. ( |22| . They 
are compared to the corresponding soft-photon limit (curve) using 
Eq. ( |23| . The data correspond to Z = 1 and to a laser probe with 
uj = 1.55 eV. 



and normalization constant Ci^o = 2Z^/^, while the pho- 
toelectron is described by a continuum wavefunction with 



asymptotic form given in Eq. (111. In the case consid- 



ered here, A = 1 while L — 0,2. Replacing in the general 
expression of the 

?L,i,o(fc; eis + f^) 



expression of the amplitude Eq. ( 54 ) , one has 
2 



Ci.oNk dr2r2sm[kr2 + ^k.L{r2)] 
Jo 

poo 

/ drir2e-^''ie*[''''>+*''-^(''>)l sin[Kr< -f <i>^,x{r<)] 
Jo 

(56) 

and splitting the integration ranges over r> and r<, one 
gets: 

2 f°° 

TL,i^o{k;eis + il) ^ -j^CifiNk J dr2r2Sm[kr2 + ^k,L{r2)] 

/ drir^e"^''! sin[Kri -I- $„.i(ri)] 
Jo 



x{e 



i[K;r2+$„, 1(7-2)] 



-hsin[Kr2 + $K,i(r2)] / drirle-^'''e'^'''"'+'^-^''^'''^^ 

J )"2 

(57) 

From the two ri— integrals present within the braces, the 
first one is by far dominant. This comes from the presence 
of the exponentially decaying wavefunction of the ground 
state, which ensures that the integration range containing 
the origin is dominant. Thus, one can neglect the second 
term and one has: 

/ dr2r2 sm[kr2 + $fc,L(r2)]e'["''^+*-^^"=)l 
Jo 

X / dnr^e-^''! sin[«;ri -^$^,i(ri)] (58) 
Jo 
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We remark that the ri-integral contained in this expression 
is real: J dri... £ R and its precise value does not affect 
the overall phase of the amplitude which becomes: 

arg[Ti_i^o(fc;eis + ^)] « 

poo 

arg{- / drrsin[fcr + $fc,L(r)]e*[''''+*'''i('')l} (59) 



or, making explicit the expressions of the Coulomb phases 
writing the sine function under its exponential form and 
keeping only the r— dependent terms in the integrand, the 
amplitude can be expressed in terms of the integrals J± 



defined in Eq. (19). Following the same line of reasoning 



as the one followed in Sec 3.2 one can neglect the contri- 
bution of the J+ integral, so that: 



arg[Ti^i,o(A:;eis 
(2fc) 



iZ/k^—i<JL(k) 



: arg{(-)^(2K)'^/''e''"i('') 



where the last integral represents a Gamma function of 
complex argument times algebraic factors, so that one re- 
covers the expression of the phase of the radial amplitude 



given in Eq. (21 1 



This expression, which has been established for an ini- 
tial Is hydrogenic state with nuclear charge Z, remains 
valid for any s bound state with an exponentially decaying 
wavefunction. As already mentioned, for non-hydrogenic 
systems, the Coulomb phase-shifts have to be replaced by 
the relevant scattering phase-shifts . 

For the sake of reference, we give the complete expres- 
sion of the asymptotic form of the amplitude, including 
real factors: 



T, 



L,1.0 



(A:;ei, +n) 



2 y K k> 



{2kyz/k 



(k- kY + tZ(l/K-l/k) 



fOO 

X / drr^e" 
'o 



sm nr 



['tr + $„,i(r)] (61) 
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